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, '^Z Abstract. We characterize the interaction between a single atom or similar 

microscopic system and a light field via the scattering ratio. For that, we first derive the 

CO ' electrical field in a strongly focused Gaussian light beam, and then consider the atomic 

, . response. Following the simple scattering model, the fraction of scattered optical power 

\^ I for a weak coherent probe field leads to unphysical scattering ratios above 1 in the 

OO ' strong focusing regime. A refined model considering interference between exciting and 

^. , scattered field into finite-sized detectors or optical fibers is presented, and compared 

\l ' to experimental extinction measurements for various focusing strengths. 
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1. Introduction 

Atom-light interaction at the single quantum level plays an important role in many 
quantum communication and computation protocols. While spontaneous emission can 
provide a natural transfer of atomic states into photonic qubits, strong interaction of 
light with an atom is needed to transfer a photonic qubit into internal atomic degrees 
of freedom as a stationary qubit. This process is essential to implement quantum light- 
matter interfaces [IIEIIS], unless post-selection techniques are used |1]. 

The common approach to achieve this strong interaction pursued for a long time 
is to use a high finesse cavity around the atom, in which the electrical field strength 
of a single photon is enhanced by multiple reflections between two highly reflective 
mirrors, resulting in a high probability of absorption. Another approach to increase 
the interaction between an atom and a single photon is simply to focus the light 
field of a single photon down to a diffraction limited area, motivated by the fact that 
the absorption cross section of an atom is on the order of the square of the optical 
wavelength. Recent theoretical research on this matter predicts that the absorption 
probability may reach the maximal value of 100% for dedicated focusing geometries [5]. 
In this paper, we study the interaction strength between a two-level system and a tightly 
focused weak coherent light Gaussian beam, which is simpler to prepare. 

Such a system has been theoretically investigated by van Enk and Kimble ^] and 
they concluded that one can expect only a weak interaction. An experiment on single 
atom absorption has been carried out a long time ago in the weak focusing regime 
[7j, but recent experimental results with single molecules [8] and atoms [9] showed an 
interaction strength which exceeded these theoretical predictions by far. In this paper 
we extend the original theoretical model such that it is applicable in the strong focusing 
regime and provide experimental data on the extinction by a single atom for various 
focusing parameters. Extrapolating from there, we find that the interaction strength in 
between light and an atom can indeed be very strong for realistic focusing geometries. 

The paper is organized as follows: In Section [2], we explain how we quantify the 
interaction strength between an atom and a weak coherent light field, and set out the 
basic problem. In Section [3l we calculate the field strength at the focus of an ideal 
lens by considering a Gaussian incident beam for the strong focusing regime. Using 
this 'ideal' focusing field developed in this section, we then obtain an expression for the 
scattering ratio in Section HI and for the extinction of a focused light beam by a two- level 
system in various geometries in Section [51 The theoretical prediction is compared with 
our experimental results in Section [6l 

2. Basic problem 

The system that we investigate is a single two-level atom localized in free space 
illuminated by a focused weak monochromatic light field (probe) with an incident power 
Pin. The interaction strength of the probe with the atom is directly related to the 
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Figure 1. the electrical field E-m of a collimated beam with Gaussian profile is 
transformed into a focusing field Ep with a spherical wave front by an ideal thin 
lens with a focal length /, leading to a field amplitude Ea at the location of an atom. 



fraction of power scattered by the atom. Therefore, it seems reasonable to quantify the 
interaction strength with the ratio of the scattered hght power P^c to the total incident 
power Pin, i.e. 

^.c:=^. (1) 

-'in 

To prepare an atom into a clean two-level system, it is convenient to use optical 
pumping with circularly polarized light - the optical transition therefore will also be 
driven by circularly polarized light. The light field itself should have a well-defined 
spatial profile before it is focused onto the atom with a lens. A circularly polarized, 
collimated Gaussian beam propagating along the z axis (see figure [1]) will therefore be 
the starting point for our work. Its electrical field strength before the lens is given by 

El 



E{p, t) = ^ [cos{ujt)x + sm{ut)y] e"^ ^""^ 



(2) 



where p is the radial distance from the lens axis, w^ the waist of the beam, x, y are the 
unit vectors in transverse directions, and E-^ is the field amplitude. The beam carries a 
total power of 
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-eoTTcE^w-^ , 



(3) 



where eo is the electric permittivity of vacuum, and c the speed of light in vacuum. Due 
to the rotational symmetry the field on the lens axis is always circularly polarized. For 
an atom that is stationary at the focal point of the lens, the electric field can thus be 
written as 



E{t) = —^ [cos{ujt)x + sin(ct; 
v2 



(4) 
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where Ep^ denotes the amphtude of the field at the focus. In the long wavelength limit 
the atom only interacts with the field at the location of the atom. For a field that is 
resonant with the atomic transition and with an intensity much below saturation, the 



power scattered by a two-level atom is [TO] (see [Appendix A| for more details) 



2Z?2 



_ 3eocA Ej^ 
^- - ^ ' l^J 

leading to a scattering ratio of 

which is exact under weak and on-resonant excitation. To evaluate the the scattering 
ratio Rsc and therefore the interaction strength, one needs to know {Ep^/ Ei)"^. 

For a weakly focused field where the paraxial approximation holds one finds that 

where Wf is the Gaussian beam waist at the focus. This leads to 

3A^ 
R^, ^ — ^ = 3u^ , (8) 

with the focusing strength 

u := Wl/I (9) 

used to fix the focal waist, Wf = X/{ttu). With a Gaussian focal spot area A = TTwf/2, 
the scattering ratio can also be expressed as Rsc — cTmaxM, where (Xmax = 3A^/27r is 
the absorption cross section of a two-level system exposed to a resonant plane wave. 
However, for strongly focused light, the paraxial approximation breaks down, and we 
need other methods to find {E^/ Ei^Y . 

3. Electrical field in a tight focus 

The paraxial approximation breaks down for strongly focused beams both in the 
expression of the electric field just behind the lens, and in the propagation of this field to 
the focus. An approach to overcome the propagation problem was reported by van Enk 
and Kimble in [6]. Their lens model, however, applies only to the weak focusing regime. 
In the following, we present a lens model to overcome this limitation and propagate the 
optical field behind the ideal lens into the focal regime and investigate the focal field 
using their technique numerically. We obtain a closed expression for the electrical field 
in the focus using the Green theorem for the propagation. 

To simplify the expressions in this section, we express the electrical field in 
dimensionless units, so the electrical field strength of the coUimated Gaussian beam 
entering the focusing lens is given by 

Fi, = e^e~P"l< , (10) 

where e+ is one of the circular polarization vectors e-t = (x ± iy)/ \/2. 
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3.1. Model of an ideal lens 

An ideal converging lens converts a beam with a plane wave front into one with a 
spherical wave front which converges towards the focal point F. Therefore, it can be 
modeled as a phase plate modifying an incoming field Fi^ with a radially dependent 
phase factor ip{p) into 

FF = <^(p)Fi„. (11) 

In paraxial optics, a convenient analytical treatment of Gaussian beams can be obtained 
assuming a parabolic phase factor, 

y.pb(p) = 6-^'"'/'^ ; (12) 

this was adopted in [6j. However, the conversion of a plane into a spherical wave front 
corresponds to a phase factor of 



^sp(p) = e-^'=v/^^^, (13) 

which is only approximated by equation (TT2l) . On top of this, multiplication of a 
incoming field with such a phase factor leads to an electrical field which is not compatible 
with Maxwell equations, since the polarization vector for p > is not tangential to the 
wave front anymore. 

In view of this, we have to change the local polarization with three requirements 
in mind pjj: (i) A rotationally symmetric lens does not alter the local azimuthal 
field component, but tilts the local radial polarization component of the incoming field 
towards the axis; (ii) The polarization at point P (see Figure [1]) after transformation 
by the lens is orthogonal to the line FP; (iii) The power flowing into and out of an 
arbitrarily small area on the thin ideal lens is the same. These requirements determine 
completely the focusing field right after the lens. With the input field in equation (fTOl) . 
one finds (see Appendix B for details) 



f^ , , ,, 1 fl + cos9^ sin^e*'^. cos^-12,va. 



X exp {-p^/wl) exp -ik^ p^ + /^ , (14) 

with d = arctan(p//). In particular, the factor l/vcos^ is needed in order to meet 
requirement (iii). 

3.2. Numerical propagation of the field to the focus 

The optical field with a converging wave front directly behind the lens needs to be 
propagated into the focal region to arrive at a field strength of the light interacting 
with the microscopic system. Various methods can be applied for this purpose. The 
one implemented in [6] projects the focusing field Fp on an orthogonal set of modes 



F^, p = (kt,s,m) with cylindrical symmetry (see Appendix C for details). This 
decomposition reads 

Ff = ^/.^F^ (15) 
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where the expansion coefficients k^ are given by 

A I. r, 1 j sk + K ( l + cos9 \ V2ktfsme\ 

n, = S^^nk, j^ dp p-j= 1^— [-^^) Uhp) + ^^ ^^J Ji( V) 



+ ^ ( 2 1 Mhp)^exp 



-tk,/^Tp - 4 



(16) 



with 6 = tan~^(p//). The Kronecker symbol 5^1 reflects conservation of angular 
momentum under the lens transformation p^ [T3] . The projection integral has no 
analytic solution, so the coefficients have to be evaluated numerically. 

The (dimensionless) field components in the three polarization components e±,z at 
any point behind the lens are superpositions of contributions from different modes, 

F+(p, ^,z)=Yl f dkt ^'-^^Mhp)e^''^^K„ (17) 



s=±l' 



F,(p,0,z) = J2 [ dh (-^)^^Ji(Me^'^^e%, (18) 

F_{p,<p,z) = Y, I dk, l-'±^J^{k,p)e^'^^e'^f>K,. (19) 

We now evaluate the field components for different regions with this method: ffist 
directly behind the lens, then on the optical axes to near the focus, and finally in the 
focusing plane near the focus. 

3.2.1. Focusing field reconstruction As a consistency check, we ffist evaluate the field 
components right after the lens (i.e., for z = — /) for a reasonably strong focusing field 
with u = 1.56 corresponding to wl = 7 mm for / = 4.5 mm. The relative difference 
between the reconstructed and original field is less than 10~^, a bound limited by our 
numerical accuracy. A linear combination of the field modes p in the form of equations 
( IT71) to ( TT9l) is compatible with the Maxwell equations, so since there is no significant 
difference between the original and reconstructed field, the choice in equation (fT^ for the 
focusing field is compatible with Maxwell equations even for strong focusing parameters. 

3.2.2. Field along the optical axis The solid line in figure [2] shows the dimensionless 
intensity l-F+P for / = 4.5 mm, A = 780 nm and wl = 1.1 mm [u = 0.244), with a clearly 
peaked distribution centered at the focus Az = and a depth of field, defined as the 
full width at half maximum (FWHM), of about 9.5 /im. This result is still very close to 
the much simpler paraxial approximation of a Gaussian beam, with a depth of field of 
2X/{nu^) =8.31/im. 

For comparison, we show the result for a focusing field using a parabolic phase 
factor (y9pb only, following [6]. The spherical aberration there displaces and spreads the 
focus, and significantly reduces the maximal intensity in the focal point F. This problem 
becomes even more serious for a larger input waist u7l- 
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Figure 2. Dimensionless intensity jF+p along tlic optical axis for the full focusing 
field model according to equation (|f4p . the paraxial Gaussian beam model, and for 
comparison with a parabolic phase factor for the lens only. The latter leads to spherical 
aberration, manifesting in a spread of the focal area and a shift towards the lens. 



3.2.3. Field in the focal plane We now examine the field near the focus in more 
details. Figure [3] shows the field on the focal plane for different focusing strengths. 
For this, we choose different input waists wl, but keep /=4.5mm and A=780nm fixed. 
For comparison, we also show the result for a focusing field according to the paraxial 
approximation. 



par ax 



!f^; , p-pV-i 



Wi 



-t^e 



(20) 



with a paraxial focal waist Wi = /A/vtwl- For weak focusing {u = 0.022, wl = 0.1 mm) 
in figure [3^, |F+| overlaps completely with the paraxial prediction with negligible 
\Fz\ and \F^\. For an initial waist wl = 0.3mm corresponding to u = 0.067 and 
Wf ^ 3.7/im (about 5A), discrepancies between the paraxial approximation and the 
extended model start to appear (Figure [3]d). With increasing wl, the z- and e_ polarized 
field components become stronger for p > 0, but an atom localized on the optical 
axis still only experiences a e+-polarized field. Figure [SJi shows the focused field that 
maximizes \F^\ for the parameters in our model. It is obtained with an incident waist 
wl = 10 mm {u = 2.22). An increase of the incident waist beyond that does not reduce 
the focal spot size any further due to the diffraction limit. Instead, more energy is 
transferred to the other polarization components, thus decreasing the magnitude of the 
F^. 
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Figure 3. Field amplitudes at the focus for different focusing strengths. All plots are 
for a focal length of 4.5 mm and wavelength of 780 nm. 



3.3. Analytical expression for the field in the focal point 

An alternative method for propagating the focusing field directly behind the lens is 
offered by the Green theorem. For given electrical and magnetic fields E[f') and B{f') 
on an arbitrary closed surface S' that encloses a point r, the electrical field at this point 
is determined by [H] 

E{r) = (f dA' {ike \n' x B{r')\ G{f,r') + \n' x E{f')\ x V'G(f,r") 



+ \n' -Eif') 



V'G(f,f')}, 



(21) 



where n' is the unit vector normal to a differential surface element dA' and points into 
the volume enclosed by S', and G{f,f') is the Green function given by 

„ik\r—r'\ 

G{r,f') = J" (22) 

47r|r — r \ 

If point r is the focus of an aplanatic focusing field, then the local field propagation 
wave vector k ' at any point f' always points towards (away from) point f for the 
incoming (outgoing) field in the far field limit, i.e. when \f— f'\ ^ A. In this limit, one 
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X E{f'), 
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VG -^ -ik'G before the focus , VG 

In the far field hmit, Equation (|2T1) reduces to 



ik 'G after the focus . 



(23) 
(24) 
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-2i dA' 
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+ 2i dA' 



n ■ k 



E{r')G{r,r') 
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Eir') 



k'G{f,f') 



(25) 



Here the surface S' is divided into two parts, where Sm is the one before the focal plane, 
and S'af is the surface after the focal plane. The second term in Equation ( l25l) vanishes 
if we choose S'af to be an infinitely large hemisphere centered at the focus, since in 
this case n' is perpendicular to E{f') at all points on S'af for an aplanatic field. If we 
choose S'bf as an infinitely large plane that coincides with the ideal lens and adopt the 
dimensionless incident field in Equation ( TT^ . we get 

-tk^ 

'0 



F(0, 0,2 = 0) 



' p(/ + v/7^+7) . p\^ 

dp — fj,o , _9N.^M — exp(— ^) e+ 



{P + P^f" 



w. 



(26) 



which has an analytical solution 



F(0,0,. = 0) = -l^^eV«^ 
4 u 



r(-Z'A) + V^r(i,^; 

u 4 u^ 4 M^ 



e+ 



(27) 



with the incomplete gamma function r{a,b) = j^ f^^^e^^ dt , and u = wi/f as in 
Equation iQ. The results obtained with the mode decomposition method agree with 
this expression within computational errors of about 0.1%. The —i reflects a Gouy 
phase of — 7r/2 fT5]. 

We now restore the field dimensions by multiplication with the amplitude in the 
center of the coUimated Gaussian beam, which can be expressed by the optical power 
according to equation ([3]), 
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resulting in an electrical field amplitude in the focus of 
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(29) 



with purely circular polarization. The focal field thus only depends on the input power, 
the optical wavelength and the focusing strength u. 
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Figure 4. Scattering ratio Rsc as a function of the focusing parameter u — wi,/ f with 
the fuU focusing field, in paraxial approximation and with a parabolic wavefront. 



4. Study of the scattering ratio 

The electrical field amplitude Ea at the focus for a given excitation power now allows 
to determine the fraction Rsc of the optical power scattered away by a two-level atom 
or similar microscopic object according to equation ([5]). With Equation (!29|) . we arrive 
at 

2 



-Rs, 



4m3 



.2/«2 



4 M^ 4 M"^ 



(30) 



In Figure H] we show this quantity as a function of the focusing strength u. A 
striking feature of this plot is that, for u large enough, Rsc exceeds the value of 1, as if 
more light would be scattered than was incident. However, Rsc cannot be interpreted 
as a scattering ratio anymore if the solid angle subtended by the excitation field is not 
negligible: in the strong focusing regime, interference between the exciting field and the 
scattered field must be taken into account. In fact, the physical bound (Bassett limit) 
is -Rsc < 2 [TTj. For our focusing model of the Gaussian beam, we predict a maximal 
value of -Rsc = 1-456 for a focusing strength u = 2.239. 

For reference, we also show -Rsc for focal fields derived under the paraxial 
approximation and for a parabolic wave front model. All models agree in the weak 
focusing regime, and -Rsc can reach values on the order of 1, which indicates that by 
strong focusing, one can accomplish a strong interaction between a light field and a 
single atom. 
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Figure 5. A transmission measurement setup with an atom at the focus of a lens. 
The transmitted power is a result of interference between the scattered light and the 
probe for coherent scattering. 



5. Extinction as a measurement of scattering 

The parameter Rgc we have just discussed can be used as a figure-of-merit for scattering 
experiments. In this section, we go beyond this one-parameter description to provide a 
detailed model for the experiment we performed p. In this experiment, we measured 
the extinction of the transmitted light beam due to scattering. 

Figure [5] illustrates a simple transmission setup with an atom located at the focus of 
two confocal lenses, where the second lens collects all the excitation power if no atom is 
present at the focus. The actual measured transmission T depends on the area covered 
by the power detector after the second lens, and can be obtained by considering the 
interference between the incident probe field and the field coherently scattered by the 
atom. For a given focusing and collection geometry, the extinction 

-Mn -'out 



1-T 



Br 



(31) 



is maximal for a weak incident probe field in resonance with the transition in the two- 
level system. The total field at any place is a superposition of the focusing field exciting 
the atom, and the scattered field: 

Et(f)=^F(r) + ^sc(r). (32) 

The spatial dependency of the scattered field Egc is that of a rotating electrical dipole, 
with an amplitude proportional to the exciting electrical field amplitude E^, and the 
total power contained in this dipole radiation must match Equation ([5]). Far away from 
the dipole (r ^ A), this scattered field takes the form 

3^^gi(fcr+7r/2) 



Esci^ 



2kr 



[e+ 



ie+ ■ r J rj , 



(33) 



where f is the radial unit vector pointing away from the scatterer [10]. The 7r/2 phase 
refiects the fact that the dipole moment of the atom lags the field Ea by tt/2 at resonance. 
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The focusing field Ep close to the lenses at z = ±/ takes the form 

El f It cos 9^ _sin6'e*'^„ , =F cos 6^-1 2i</,. 



^^^'■'''-^^'^-Tm^^'^^^r'^^^'''- 



X 



exp(-pVwL) exp ±i{k^p^ + P- 7r/2) , (34) 



where 9 G [0, vr] is the polar angle between the —z direction and a point (p, 0, z) as in 
figure m The phase is adjusted such that the electric field amplitude Ey^ at the focus is 
real. 

The excitation field and the forward scattered field interfere destructively, as was 
shown first for the case of an incident plane wave [161 UHl EH], and more recently for 
arbitrary incident fields [20] with the help of vectorial multipole expansions 



5.1. Energy flux through transverse planes 

The optical power Pout arriving at a detector behind the atom or a collection lens can 
be evaluated from the superposition of fields via the time averaged energy flux through 
a plane S with a flxed z, 

Ps = ^ halEtX B*] ■ dl (35) 



^ Js ^ ■' 

where dA is a differential area element of the surface S and 3fJ(x) denotes the real part 
of X. 

Far away from the focus, the electromagnetic fleld can be locally approximated 
by a plane wave such that B = k x E/c, where k = k^c, kp is a dimensionless unit 
vector parallel to the local fleld propagation direction. Both Ep and E^c have spherical 
wave fronts, i.e., the local propagation directions are parallel. Before the focus we have 
ksc = —kp, while after the focus we have ksc = kp. With these fleld properties, and with 
the local transversality, kp ■ Ep = and kp ■ E^c = 0, the power through the two planes 
can be expressed with electrical flelds only, 

P.=±f = ^ / ^{Ep-E^±E,,- 4; + ise ■E*p±Ep- 4;} kp ■ ma . (36) 

The two flrst terms represent (i) the power of the excitation fleld, (ii) the power of the 
scattered fleld, while the third and fourth term represent (iii) the interference term. 
The contribution (i) to Pz=±f is simply the input power, 

^ y ^[Ep-WpYkp- zdA = ^eoTTcElwl = P,^ , (37) 

and the contribution (ii) in these planes, 

\P.^±f,s.\:=fl^^^{E..-K}k.-MA=^-^^^^ (38) 

where Psc is the scattered power as deflned previously in Equation ([5]). 



=±f,m 

/z=±f 
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The interference contribution (iii) vanishes for z 
is purely imaginary, whereas for 2; = +/ we get 



/ because ( E^c ■ Ep - Ep ■ E* 



P.-. 



=+f , int 



=±f 



E; + Ef- E*^ \ kp ■ zdA 



SneocE^ELvl r p{f + ^/f^T?) 



2k 



(/2 + p2)5/4 



exp( 



^)dp. 



w 



(39) 



(40) 



The negative sign, which conies from both the Gouy phase in the incident field (Equation 
( 1M|) ) and the phase difference between the dipole and local field (Equation ([2SD), reveals 
that the scattered light and the incident light interfere destructively after the focus |20j . 
This integral can be solved in the same way as Equation fl26|) . leading to 



P. 



z=+f , int 



^eocX^El 



-Ps. 



Thus, the power flowing through both planes z = ±/ is the same, 

Pz=±f 



p 

p. — 111 



(41) 



(42) 



This indicates what we mentioned above: a value Rgc > 1 does not violate energy 
conservation, the physical bound being rather Rsc < 2. 

For later on, we should define a measurable extinction as the difference of the 
transmitted power with and without the atom, divided by the power transmitted without 
the atom at a location behind the atom, e.g. at 2; = +/: 

_ Pz=+i, in — Pz=+{ 
Pz=+i, in 

-* z=+f , int -^ z=+f,sc 

Pz=+i, in 

If the collection lens at z = +/ is infinitely large, it takes a value of 

-* SC / ^ -Kgc 

^ ~ ~P ~ ~2~ ' 

5.2. Extinction observed with a detector /lenses with finite diameter 



(43) 



(44) 



Realistic lenses will have a finite size, thus only partly transmit the excitation and 
scattered light. We now estimate how this obstruction affects the relation between an 
observed extinction and the inferred scattering ratio -Rsc- 

The first effect of a finite lens aperture radius po of the first lens is a reduction 
of the field at the focus. The Green theorem method for evaluating the focal field via 
equation fl26l) still has a closed solution for a finite radius. 
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(45) 
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with u = wl/ f as in Equation (^, and similarly v := po// as half the f-number of the 
lens which is related to the numerical aperture NA of the lens via NA^ = f^/(l + v^). 
This obstruction reduces the scattered power Psc- For realistic lens sizes, however, this 
is a very small effect. For instance, for po = 2wl we find P^° ~ 0.97Pg^. 

Without the atom, the transmitted power after the collection lens is given by 



pPo _ p 

-'^f=+z,m ~ -"^ir 



exp(- 



2p, 



2 



wl 



(46) 



The transmitted power in absence of an atom is thus very close to Pin for po > 2wl- 
The contribution (iii) of the interference terms to the forward power is now given by 

P.%f,int = -Ps'c\ (47) 

and the contribution (ii) of the scattered light by 

1 — a 

Y 



P.%f,sc =^^-^s? (48) 



with a = i±^ = (l-^)^r3^. (49) 

(1 + ^2)3/2 ^ 4 ^ 

The extinction measurable with a finite aperture lens thus is 
1 + Q, ppo I 

If the lenses fully accommodate the Gaussian incident beam, say po > 2wl, then this 

can very well be approximated by 

1 + a 
e = ^-RTc, (51) 

where P^« = PfcV^m- 

We can also quantitatively evaluate the reflectivity of a single atom in this strong 
focusing regime, which was recently found to be possibly very large [20|. The scattered 
power recollected by the input lens is also given by Equation ( HSj) . thus the 'single atom 
reflectivity' with a Gaussian beam profile is 

R = ^-R's!- (52) 

We conclude that the measured extinction presents a lower bound to the scattering 
ratio if the collection lens fully collects the probe after the focus. For a small numerical 
aperture of the collection lens, e ~ Psc as expected, whereas for a large collection 
aperture (corresponding to a small loss factor a), a reduced extinction e — * Psc/2 should 
be observed. 

5.3. Extinction observed with a detector behind a single mode fiber 

The symmetrical arrangement of the focusing and collection lens, and the typical 
preparation of a Gaussian excitation beam by an optical fiber suggests that the light 
could also be collected by a single mode optical fiber. The confinement of the light field 
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into waveguides with well-defined mode functions makes the atom-focusing arrangement 
an independent building block for 'processing' electromagnetic fields. 

The amplitude ttc of the light field collected into the optical fiber is again given 
by the sum of the excitation and scattered field, picked up by the optical fiber. This 

— * 

amplitude Oc can be obtained by projecting the field Et onto the field mode g of the 
optical fiber. This projection can be carried out with the scalar product 



^ z? \ 'J 

9,Et) := ^- 



eoc 
2 



x€S 



rix) 



■ Et{x) \ {kg-h)dA. 



(53) 



where the integration plane S is chosen such that both g and E are far away from a 
focus, kj^ is the local propagation direction of the mode function g, and n the normal 
vector on the plane S. This integration can be carried out at any convenient location as 
long as it captures the mode function. The scalar product in Equation fl53|) is written 
such that it resembles the form of the power integral in planes z = ±/ in Equation ( l36l) . 
so we can conveniently use the integrations carried out earlier. Thus, the integration 
plane S is chosen at z = +/, directly before the collection lens. 

In the experiment, the excitation mode is matched to the collecting single mode 
fiber. Correspondingly, we define the target mode function g{x, y, z) to be the same as 
that of the excitation mode of E-p in equation ( IMj) . With the normalization condition 
{g,g) = 1 we simply can set 

g = Ep/^/K. (54) 

With this normalization, the square of the projection coefficient Cc has the dimension 
of a power. Thus, the optical power of the field coupled into the fiber with a scattering 
atom present is given by 



P. 



out 



9,E, 
9,Ef 
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+ {g,E,, 



(55) 

The first scalar product is determined by the mode normalization. The second 
one, (g,Esc\ represents the projection of the scattered field onto the collection 
mode. Modulo the normalization constant a/HT, it is identical to half the interference 



contribution in Equation fl39l) . whose explicit expression was given in Equation (HI 
Therefore we find 
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(56) 



In the weak focusing regime where Psc <^ 1, this translates again into an extinction 
e ~ Rsc ■ For a focusing parameter u = 2.239, we get a maximal extinction of 
^max = 0.926. For the light scattered back into the excitation mode, we do not have 

2 

?2 



to consider the field E-p, and arrive similarly at Pback 
reflectivity of 

R = RlJA . 



9,E, 
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PinRi/4:, or a 
(57) 
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Figure 6. Experimental setup for measuring the extinction of a light beam by a 
single atom. AL: aspheric lens (/ ~ 4.5mm, full NA ~ 0.55), P: polarizer, DM: 
dichroic mirror, BS: beam splitter with 99% reflectivity, A/4, A/2: quarter and half 
wave plates, Fl: filters for blocking the 980 nm FORT light, F2: interference filter 
centered at 780 nm, Dl and D2: Si-avalanche photodiodes. Four more laser beams 
forming the MOT lie in an orthogonal plane and are not shown. 



6. Experiment 

In this section, we consider the results of our experiment where we measured the 
extinction of a Gaussian beam by a single ^^Rb atom with different focusing strengths, 
and compare the results to the above theoretical model. A detailed description of 
the experimental setup is reported in |9j and shown in Figure |6l Two aspheric lenses 
(/ = 4.5 mm) are mounted in a UHV chamber in a confocal arrangement. A single ^^Rb 
atom is localized in a far-off resonant dipole trap (FORT) that is formed by 980 nm 
light at the focus of the lens pair. A probe beam is delivered from a single-mode fiber 
and focused onto the atom by one lens, and picked up by the other one. The confocal 
arrangement ensures that all of the incident probe power is collected in the absence of 
an atom, thus implementing the scheme discussed in the previous section. We use a 
circularly polarized probe to optically pump the atom into a closed cycling transition. 
After allowing some time for optical pumping, we measure the transmission of the 
probe beam that is defined as the ratio of count rates at detector Dl when the atom is 
present in the trap, to the count rate when the atom is absent. Such a measurement 
is carried out for different probe frequencies to obtain the transmission spectrum of a 
single Rb atom. The spectrum is fitted to a Lorentzian with the resonant frequency, 
the full width at half maximum (FWHM) of the spectrum and its minimum value Ty^m 
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on resonance as parameters. We obtained spectra for four different input waists of the 
probe, thus measuring extinctions for different focusing strengths. The observed FWHM 
never exceed 7.7 MHz, which is close to the natural linewidth of the optical transition 
(6 MHz), so we conclude that the atom was successfully kept in a two-level system. 
For each probe frequency, the probe power was adjusted such that the atom scatters ~ 
2500 photons per second, way below saturation. The properties of various transmission 
spectra obtained with different probe incident waists are summarized in Table [TJ 



Table 1 . Summary of transmission spectra of the probe for different focusing strengths 
u. Listed are wi^: incident waist of the probe; W{ and wb: estimated paraxial focal 
waists of the probe beam and FORT, respectively; e and W: maximum observed 
extinction value and FWHM of the transmission spectrum; Rsc- scattering ratio for 
the focusing parameter; ethco: expected extinction according to Equation (|56|) . 



Wh (mm) 


u 


Wf(^m) 


WB (A^m) 


emax (%) 


W (MHz) 


Rsc 


ethco(%) 


0.5 


0.11 


2.23 


2.0 


2.38 ± 0.03 


7.1 ± 0.2 


0.0362 


3.58 


1.1 


0.24 


1.01 


2.0 


7.2 ±0.1 


7.4 ± 0.2 


0.1606 


15.41 


1.3 


0.29 


0.86 


1.4 


9.8 ±0.2 


7.5 ± 0.2 


0.2157 


20.40 


1.4 


0.31 


0.80 


1.4 


10.4 ± 0.1 


7.7 ± 0.2 


0.2449 


22.99 



We also carefully characterized the losses of the probe beam in its optical path to 
ensure that our measured extinctions are not exaggerated by interference artefacts that 
can happen when certain elements in the transmission path preferentially filter more 
probe than the scattered light [S]. From point A to point B in Figure Owe measured 
53-60% transmission without the atom in the trap. The losses are mostly determined by 
21% loss through the uncoated UHV chamber walls and 17-24% loss due to the coupling 
into the single-mode fiber at the transmission measurement channel. The coupling loss 
into the fiber increases as the input waist of the probe beam i^l increases. Almost all 
losses can be ascribed to reflections at optical surfaces, except for a 9-16% re-coupling 
loss into a single mode fiber due to mode mismatch. We are thus reasonably confident 
that our measurement is free from artefacts that may arise due to incomplete collection 
of the probe. 

In Figure [7] we compare the extinctions obtained from the experiment with values 
predicted by Equation (l56l) . Since the probe is recoupled into a single mode fiber for 
every experimental point using a different lens that matches the probe waist, and the 
coupling lens has a NA=0.55 (corresponding to w = 0.66), we can safely neglect any 
clipping, and use Equation fl56|) for estimating the extinction. 

Obviously the measured extinctions are smaller than those predicted, especially for 
a larger focusing strength. We see a few possible reasons for this discrepancy. Firstly, 
the lenses we used in the experiment may not be suffciently close to an ideal lens, since 
they were designed for a situation with an additional window in the focusing part, which 
we did not have in our experiment. Secondly, the interaction strength is significantly 
affected by the motion of the atom in the dipole trap. While we do not have an 
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Figure 7. Experimentally measured extinction (black symbols) for several focusing 
strengths u and predicted values for extinction e and Reflectivity i? for a coupling into 
single mode optical fibers. The inset shows the prediction for much stronger focusing 
parameters. 



independent measure of the position fluctuation of tlie atom in the trap, measurements 
in a similar trap showed a temperature of the atom in the trap of around 100 /iK due 
to the overlapping MOT [21], slightly below the Doppler temperature of 143 /xK of 



Rubidium. With our trap frequencies of Vp ^ 70 kHz in transverse and v^ ^ 20 kHz 
in longitudinal direction, this results in position uncertainties of Op ^ 220 nm and 
(Xz ~ 780 nm, respectively. The scattering ratio -Rsc gets reduced due to the presence of 
the atom in regions with a reduced excitation field, and due to the spatially dependent 
detuning in the optical dipole trap. In paraxial approximation, we find an approximate 
reduction of the scattering rate due to the lower average field of 

K^ ^ i?sc (1 - 2al/w^Y (1 - a'^X'/iAf)) , (58) 

which results in a reduction of 2% for wl = 0.5 mm to 23% for wl = 1.4 mm. 
The reduction in Rgc and therefore in e in our regime is about proportional to the 
temperature, so a doubling of the temperature alone would explain the discrepancy 
between theory and experiment already. On the other hand, the contribution due to a 
spatial variation of the resonance frequency is less than 1%. Additionally, the presence 
of the atom away from the focal point reduces the efficiency of the optical pumping. 
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However, the observed extinction ratios still exceed the prediction by the parabolic wave 
front model in [6]. 

The inset of Figure [7] extrapolates the extinction we could expect for much stronger 
focusing; as mentioned earlier, the maximal extinction should reach 92.6%; whether such 
lenses can be manufactured, however, remains an open question. We also depicted the 
reflectivity of the fiber-atom-fiber system, which should reach values strong enough 
to be detected in an experimental setup. Extinction and reflectivity don't match in 
this configuration, which means that there is still a significant amount of light which is 
neither transmitted nor reflected back into the optical fiber. Such losses are unavoidable 
for the described coupling scheme, which still places the fiber-atom-fiber system in 
disadvantage to a atom-cavity system in terms of success probability of scattering into 
known modes. If it would be possible to achieve a better overlap of photonic modes 
with the dipole transition, such losses should be reduced. 

7. Conclusion 

We have demonstrated both theoretically and experimentally that a substantial coupling 
efficiency of a light beam to a single atom can be achieved by focusing a light beam 
with a lens. By modifying the model given in [6], we have constructed a focusing 
field compatible with Maxwell equations that is suitable for the strong focusing regime. 
High values for the extinction of light (up to 92%) by a two-level atom stationary at the 
focus under the assumptions of weak on-resonant coherent probe are predicted. Within 
the limitations of our current trap, our experimental results confirm the possibility of 
observing a substantial extinction already for relatively weak focusing. 

The measured extinction depends on the particular collection configuration of an 
experimental setup. It is thus not a fixed quantity for a given incident field. As such, 
the scattering ratio as defined by Equation ([6]) is a better quantity to characterize the 
interaction strength between a weak coherent field and an atom in free space, even 
though it loses a simple physical interpretation in the strong focusing regime. These 
results may also be of interest for experiments with single molecules P, [25] and quantum 
dots 1261. 
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Appendix A. Power scattered by the atom in a coherent hght field 



We briefly repeat the resuhs for hght scattered by a two-level atom following |10j . 
The steady-state population P22 of the excited state in a two-level atom exposed to a 
monochromatic field can be obtained from the optical Bloch equations: 

^'' ~ 52 + 1^172 + r2/4 ^ '^^ 

Therein, T is the radiative decay rate of the excited state, 

and Q = -E'aMi2|/^ is the Rabi frequency for a vanishing detuning S = uj — uji2 of the 
driving field with respect to the atomic transition frequency uju- Therein, \di2\ is the 
electrical dipole moment of the atom. 

The optical power scattered by this atom is simply the product of energy splitting, 
decay rate and population of the excited state: 

^sc = p22^f>^12 (A. 3) 

For weak {Q ^ F) on resonant {S = 0) excitation, the scattered power becomes 

Psc = ^^^. (A.4) 

for an excitation field amplitude Ea at the location of the atom. 

Appendix B. Transformation of local polarization by the lens 

To obtain the local polarization of the focusing field in Equation (fT4|) . we consider a 
point P(p, (f), z) before the lens and an incident light field with polarization 

X + iy 

or in the cylindrical basis, 

72" ' V2 

where p = cos (p x + sin (p y and (p = — sin (p x + cos (p y are two unit vectors along 
the radial and azimuthal directions respectively. The ideal lens leaves the azimuthal 
component unchanged but tilts the radial component such that the local polarization of 
the field right after the lens is perpendicular to the line FP in Figure [1] (F is the focus 
point), that is: 

/cos6'e*'^ sin6'e*<^ \ ie''^ - 

1 -I- cos 9 „ sin 6'e*'^ ^ cos 9 — 1 ^ • , „ 
2 + V2 2 

where 9 = arctan(p//) and e_ = (x — iy)/\/2. 



(B.l) 



e^ ie 

ei„ = -=p+— ^0, (B.2) 
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Appendix C. Decomposition of a field into modes with cylindrical 
symmetry 

For completeness, directly following [6], we outline the main properties of the cylindrical 
modes F^, which form a complete orthogonal set to compose an electric field that satisfies 
the source-free Maxwell equations, 



E{t) = 23? 



5^a.F,e^-* 



(C.l) 



where the summation over z/ is a short-hand notation for 

Y.--= j dk j dk^Y^Y.^ (C.2) 

V s m 

and a^y are arbitrary complex amplitudes. The modes are characterized by four indices 
u := (/c, kz, m, s), where A; = ^ is the wave vector modulus, k^ = k ■ z the wave vector 
component in z-direction, m an integer- valued angular momentum index, and s = ±1 
the helicity The dimensionless mode functions Fj, in cylindrical coordinates (p, z, 0) 
given in [27] are 

F^(p, z, 0) = — -G{k, kz, m + l)e_ + — -G{k, k^, m - l)e+ 

47r fc 4:71 K 

-i- — -^G{k,kz,m)z, (C.3) 

4tt k 

where kt = \/k'^ — k1 is the transverse part of the wave vector, e± = (x ± iy)/ \/2 are 
the two circular polarization vectors, and 

G(fc, A;,, m) = JUktp)e"'''e'"'t (C.4) 

with Jm the rri-th order Bessel function. As we are interested in a monochromatic beam 
with a fixed value oi k = 27r/A propagating in the positive z direction {k^ > 0), the 
set of mode indices is reduced to p := [kt, m, s) where, for convenience, kt is taken as a 
mode index instead of k^. Now, we introduce the notation 

fi s m 

— * 

for a complete summation over all possible modes. For a fixed k the modes F^ are 
orthogonal with respect to the scalar product 

F„ ^,')= [ F;{f) ■ F,if) dS = 5{h - K)6mm'^,s'/i27rkt) , (C.6) 

where S* is a plane perpendicular to the z axis. This scalar product can thus be used to 
find the amplitudes of the modes p, fi' in an arbitrary electric field compatible with the 
Maxwell equations. 
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